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(T) ' ABSTRACT 

■ In the context of the study of the Integrated Sachs Wolfe effect (ISW), we construct a 
^ | template of the projected density distribution up to redshift z — 0.7 by using the Luminous 

. Galaxies (LGs) from the eighth Sloan Digital Sky Survey (SDSS) data release (DR8). We use 

^-—l ■ a photometric redshift catalogue trained with more than a hundred thousand galaxies from the 

| Baryon Oscillation Spectroscopic Survey (BOSS) in the SDSS DR8 imaging area covering 

. ■ nearly one quarter of the sky. We consider two different LG samples whose selection matches 

that of SDSS-III/BOSS: the low redshift sample (LOWZ, z e [0.15,0.5]) and the constant 
j_j , mass sample (CMASS, z e [0.4, 0.7]). When building the galaxy angular density templates we 

use the information from star density, survey footprint, seeing conditions, sky emission, dust 
extinction and airmass to explore the impact of these artifacts on each of the two LG samples. 
In agreement with previous studies, we find that the CMASS sample is particularly sensitive to 
Galactic stars, which dominate the contribution to the auto-angular power spectrum below ( = 
7. Other potential systematics affect mostly the very low multipole range {£ e [2, 7]), but leave 
fluctuations on smaller scales practically unchanged. The resulting angular power spectra in 
the multipole range I e [2, 100] for the LOWZ, CMASS and LOWZ+CMASS samples are 
compatible with linear ACDM expectations and constant bias values of b = 1.98 ± 0.11, 
2.08 ± 0.14 and 1.88 ± 0.11, respectively, with no traces of non-Gaussianity signatures, i.e., 
/^° L cal = 59 ± 75 at 95 % confidence level for the full LOWZ+CMASS sample in the multipole 
range £ € [4, 100]. After cross-correlating WMAP-9yr data with the LOWZ+CMASS LG 
projected density field, the ISW signal is detected at the level of 1.62-1.69 <x. While this 
result is in close agreement with theoretical expectations and predictions from realistic Monte 
Carlo simulations in the concordance ACDM model, it cannot rule out by itself an Einstein- 
de Sitter scenario, and has a moderately low signal compared to previous studies conducted 
on subsets of this LG sample. We discuss possible reasons for this apparent discrepancy, and 
point to uncertainties in the galaxy survey systematics as most likely sources of confusion. 

Key words: cosmology: observations - cosmic microwave background - large-scale structure 
of the Universe 
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1 INTRODUCTION 

The growth and evolution of density and potential perturba- 
tions in an almost perfectly isotropic and homogeneous uni- 
verse can be described in General Relativity by the fam- 
ily of Friedmann-Robertson- Walker (FRW) metr ics (see, e.g., 
iLifshitz & Khalatnikodl963l ; lMukhanov et alj 19921) . Despite their 
requirements on isotropy and homogeneity, these models leave 
room for the study of small clumps or inhomogeneities that eventu- 
ally give rise to the complex structure of today's universe. Indeed, 
under the assumption that initial fluctuations in the metric are small, 
it is possible to write a linear theory of perturbations describing the 
interplay and ev olution of small inhomoge neities in a FRW back- 
ground (see, e.g. jMa & BertschingeJl995l and references therein). 

In particular, FRW models show that, on the large scales 
within the horizon where the Poisson equation holds, the growth 
rate of the gravitational potentials (D$(t)) is proportional to the 
time derivative of the ratio D p (t)/a(t), where a(t) is a FRW met- 
ric scale factor describing the expansion of scales in the universe, 
and D p {t) is a function describing the (linear) growth of matter den- 
sity in-homogeneities. In the special case of an Einstein-de Sitter 
universe D p (t) = a(t) and hence gravit ational potentials on large 
scales should remain constant (see, e.g. jMa & Bertschinger|[l99a . 
and references therein). However, deviations from this scenario will 
make D p (t) ± a(t) and thus variations in the large scale gravi- 
tational potentials will arise. Photons of the Cosmic Microwave 
Background (CMB) crossing those potentials will leave a differ- 
ent potential than they entered, and the difference will appear in the 
form of a gravitational blue/redshift that is known as the late In te- 
grated Sachs-Wolfe effect, (hereafter ISW. ISachs & Wolf J 1967b . 

If Dark Energy is responsible for the acceleration in the ex- 
pansion rate of the universe, then gravitational potential wells will 
shrink with time, and CMB photons will leave potentials that have 
become shallower during their fly-by. This results in more en- 
ergetic CMB photons emerging out of overdensities, and hence 
the ISW temperature anisotropies of the CMB should be posi- 
tively correlated to the dist ribution of gravitational potential wells. 
ICrittenden & Turokl ( 1 19961) first noticed the fact that those wells 
are probed by halos and galaxies, and that the distribution of these 
structures can be cross-correlated to the CMB anisotropies in or- 
der to distinguish the ISW component from the intrinsic CMB 
anisotropies generated at the surface of last scattering, at z — 1 , 050. 

As soon as the first data sets from Wilkinson Microwave 
Anisotropy Probe (hereafter WMAP) were released, several 
works claimed detections of the ISW at various levels of sig- 
nificance (ranging from 2 a to > 4cr) after cr oss-correlating 
WMAP data with different ga laxy surveys, (Scranton et al.l 



20031; iFosalba et alj 120031: iBoughn & Crittenden! 120041 ; 



Fosalba & Gaztanagal2004l;lNolta et al J2004I; Afshordi et alj200ll: 

r» _i i_i _ i I 1^ t\{\ c\ _ i_ . _ l I 1 /■! ■ " _ \ l 



Padmanab han et alj 120051 : ICabre et al.1 12006|: | Giannantonio et alj 
20061 ; IVielya et all l2006t iMcEwen et alj l2007h. Ho wever, other 
authors ( Hernandez-Monteagu do et al] 120061; iRassat et alj 



20071; iBielbv et alj 1201^: 



Hernandez-M onteagudcl |2010| ; 



Lopez-Corredoira et al. 
Francis & Peacock! 



2010; 



2010 



Sawangwit et alj 120101) . reported lower (and/or lack of) sta- 



tistical significance and/or presence of point source contamination 
in the cross-correlation between different galax y surv eys and 
WMAP CMB maps. In iHernandez-Monteagudol 120081) it was 
shown that the cross-correlation of the ISW with galaxy templates 
should be contained on the largest angular scales (€ < 60-80), 
with little dependence on the galaxy redshift distribution, and that 
most of such signal should arise in the redshift range z e [0.2, 1.0]. 



However, most of the surveys used in ISW-galaxy cross correlation 
studies were either shallow and/or had more anisotropy power 
than predicted by the models precisely on the large angular scales 
sensitive to the ISW. 

The issue of this power excess present at large scales, when 
computing the auto-correlation function or auto power spectrum 
in galaxy surveys, has been relevant for some time (see, e.g. , 



in galaxy survey s, has bee n relevant tor some tim e (see, e.g. , 
Ho et al. 2008; Hernan dez-Monteagudoll20ld ; lThomas et alj|201ll; 



Giannantonio et all 20121) . The shape of this excess is similar to 
that expected in models with high levels of local Non-Gaussianity 
(NG), an d thus it has motiv ated claims of NG detection in the past, 
(see, e.g-. lxia et~a l. 2010a b). Recently it has been recognized that 
the presence of systematics in the galaxy survey data m ay be largely 
responsible for this extra power (Huterer et al. 2 012b . a nd thus lat- 
est co nstraints on NG are consequently weakened, (see iRoss et al] 
l2012h . Unfortunately, the impact that these systematic effects - giv- 
ing rise to extra power on large scales - may have on the ISW 
signal has not yet been properly addressed. The last studi es on 
the subject dGiannantonio et al]|2012l ; Ischiavon et alj 2012) have 
pointed out the presence of this excess, but they did not quanti- 
fied its impact on the expected ISW signal-to-noise ratio, or galaxy 
cross-correlation analysis, i.e. either if due to NG or to systematics. 
Since this unexpected high level of anisotropy on the large scales 
is present in different galaxy surveys, i t should also affect those 
works combining different surveys (e.g. jGia nnantonio et al.ll2008t 
Ho et al. 2008; Giannantonio et al. 2012). The parallel analysis of 



Giannantonio et alT 20131) has nearly simultaneously obtained s 



ilar results in terms of ISW significance using a similar data set to 
our own. 

The Sloan Digital Sky Survey-Ill (SDSS-III. lEisenstein et al] 
201 lh Barvonic Aco ustic Oscillation Spectroscopic Survey (BOSS, 



Dawson et al] l2013h has calibrated a number of potential sys- 



tematics (star density, seeing, sky emission, etc) that have to be 
used when building templates of galaxy angular number density, 
teossetalj|20"l2l ; IHo et alj|2012h . In this work we make use of 
all the information related to potential systematics to produce cor- 
rected templates of the number density of Luminous Galaxies from 
two different BOSS samples: the low redshift sample (LOWZ, z £ 
[0.1,0.5]) and the constant mass sample, (CMASS, z e [0.4,0.7]). 
These templates are used as tracers of moderate-redshift gravita- 
tional potential wells, in order to be correlated with the CMB maps 
produced by WMAP. In Section |2]we describe the SDSS Data Re- 
lease 8 (DR8) and SDSS-III/BOSS data used for the LG template 
construction, and in Section [3] the CMB data from WMAP used 
in the cross-correlation analysis. In Section |4]we describe the the- 
oretical predictions for the level of LG - CMB cross-correlation 
to be induced by the ISW effect. The statistical methods used in 
our cross-correlation analysis are described in Section |5] and its 
results are presented in Section [6] These results are discussed in 
Section [7] where conclusions are also presented. The detailed ex- 
planation of the systematic corrections are provided in Appendix 
lAl These corrections must be taken into account when computing 
and interpreting the angular power spectrum of the LG samples, 
as described in Appendix [B] This appendix also includes tests for 
remaining residuals systematic and an assessment on the compat- 
ibility of the data to the Gaussian hypothesis (i.e., limits on the 
non-Gaussian parameter /n"™ 1 )- 

Unless stated otherwise, we employ as a reference a 
flat ACDM cosmo logical model consistent with WMAP9 
fainshawet~ai]|2012h . with density parameters O a = 0.04628, 
ficdm = 0.24022, fl A = 0.7135, reduced Hubble constant h = 
0.6932, scalar spectral index «s = 0.9608, optical depth to the Last 
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Scattering Surface r T = 0.081, and rms of relative matter fluctua- 
tions in spheres of 8 hT 1 Mpc radius erg = 0.82. 



2 BOSS LG SAMPLE 



The galaxy sample we use is based on a phot ometric se l ection 
of galaxies from the SDSS Data Release 8 jYork et al.l I2OO0I : 



lAihara et al.l 



photometry 



201 ll DR8). Th i s surv ey obtained wide-field CCD 
Gunn et al 1 ll998l l2006h in five passbands (u, g, r, i, 



z; e.g. lFukugitaetal.lll996h . The DR8 imaging data covers 14,555 
deg 2 , more than one third of which is in the Southern Galactic cap. 
This large sky area provides an unprecedented signal-to-noise ratio 
for ISW studies as compared to previous analysis in the literature, 
as shown in Section 4. 

We s elect galaxies using the same color and magnitude cuts 
defined in lEisenstein et al] bOllh . Eqs. 1-5. This includes a se- 
lection for a lower redshift (z i 0.4) sample denoted 'LOWZ' 
and higher redshift (0.4 i z £. 0.7) composed of galaxies that are 
approximately stellar mass limited and denoted 'CMASS'. These 
photometric selections yield a total of 2.4 million galaxies (after 
removing any duplicate galaxies that appear in both selections), 
one third of which are in the LOWZ sample. Photometric red- 
shifts, galaxy/star probabilities, and recommendations for measur- 
ing the_angjilar clustering for the CMASS sample were described 
by iRoss et all J201 lb and t he angular clust ering of the CMASS 
sample was studied further in lHo et aL For the LOWZ sam- 

ple we also generate phot ometric re dshifts TBoltonetal.1120121) . as 
we did for CMASS in IRoss et alj |201lh. using BO SS spectro- 
scopic redshifts as a training sample (ISmee et al . 2012). These pho- 
tometric redshifts are trained using the neural network code ANNz 
dFirfh et al.l2003l) . We remove LGs with z p i, a r < 0. 15, as such galax- 
ies are expected to have minimal contribution to the ISW signal 
whi le they can be more easily confused with other types of galax- 
ies jPareiko et al.ll2013l) . This cut removes 1.2 x 10 5 galaxies. 

The SDSS DR8 imaging area deemed s uitable for galaxy clus- 
tering studies was defined by |HoetalJ ( l2012h - areas are removed 
due to a number of reasons, s uch as image qua lity and Galactic 
extinction - using a HEALPi)fl ( Gorski et al. 2005) mask with res- 
olution Nside=1024. For the characterization analysis of the com- 
bined LOWZ and CMASS samples and the study of their auto- 
power spectrum, we propagate this mask to Nside=128 and then 
remove all data in HEALpix pixels with weight less that 0.85. This 
procedure leaves a sample of 1.6 x 10 6 galaxies occupying a 9, 255 
deg 2 sky footprint. The sample's redshift distribution, determined 
empirically from the BOSS spectroscopic redshift distribution, is 
displayed in the left panel of Fig. [I] For the cross-correlation anal- 
ysis with CMB WMAP9 data, the mask and the CMB and galaxy 
maps were further downgraded to a resolution parameter Nside=64, 
which corresponds to a pixel area of slightly below one square 
degree. In this case, in order to avoid mask border effects, only 
masked pixels with amplitudes above 0.9 were considered. In com- 
bination with the WMAP9 sky mask, this cut leaves 1.4x 10 6 galax- 
ies on 7,660 deg 2 , which represents an effective sky fraction of 
0.19. 



3 WMAP CMB DATA 

The Wilkinson Microwave Anisotropy Probe (WMAI0) scanned 
the microwave sky from 2001 until 2010 in five different frequen- 
cies, ranging from 23 GHz up to 94 GHz. The angular resolution 
in each band increases with the band central frequency, but it re- 
mains better than one degree in all bands (for the lowest frequency 
channel under consideration in our study, centred at 41 GHz, the 
square root of the beam solid angle equals 0.51°). The ISW and its 
cross-correlation to tracers of the density field arise on larger scales 
(6 > 3° or ( < 80), hence angular resolution will not be an issue in 
this context. 

The WMAP measurements have provided all-sky maps in 
which, beyond the CMB radiation, other diffuse components and 
point sources have been identified and characterised, (see, e.g., 
iBennett et alj|2012l for the latest results from the 9 years of ob- 
servations). For the data release corresponding to the first 7 years 
of observation, the si gnal-to-noise rati o was greater than one for 
multipoles E < 919 jjarosik et alj|2010 ). At the large scales of in- 
terest for ISW studies, the Galactic and extragalactic foreground 
residuals (mostly generated by synchrotron radiation, free-free and 
dust) remained b elow the 15 /uK level outside the masked regions 
dGold et alj|201 l] ). This analysis has been improved in the recent 9 
year data release ( IBennett et al.l2012h . in which an all-sky CMB In- 
ternal Linear Combination map, obtained from all bands, has been 
provided together with its covariance matrix. This analysis con- 
stitutes an optimal handling on the uncertainties imprinted on the 
CMB map by the presence of all those contaminants. 

In this context, we shall concentrate our analyses on the fore- 
ground cleaned maps corresponding to the Q (41 GHz), V (61 GHz) 
and W (94 GHz) bands, after applying a foreground mask KQ85y9 
that excludes ~ 15% of the sky. At the scales of interest, instru- 
mental noise lies well below cosmic variance and/or foreground 
residuals, and hence will not be considered any further. The ISW 
is a thermal signal whose signature should not depend upon fre- 
quency and hence should remain constant in all three channels un- 
der consideration. All WMAP related data were downloaded from 
the LAMBDA site{J 



4 THEORETICAL PREDICTIONS 

On scales large enough where linear theory applies, the Poisson 
equation relates perturbations in the gravitational potential 0^ at 
any epoch with those of the matter density contrast (6 k ) at present 



(f> k (a) = -47rGp 4 ,o 



a k 2 



(1) 



where a denotes the cosmological scale factor, D(a) is the den- 
sity linear growth factor, G is Newton's gravitational constant, k 
is the comoving wavevector and p b0 is the present background 
matter density. From this equation it is easy to see that for an 
Einstein-de-Sitter (hereafter EdS) universe, for which D(a) = a, 
the linear gravitational potentials remain constant. However, in 
all other scenarios where D(a) + a (suchs as an open universe, 
or ACDM, or any other sort of Dark Energy model) the po- 
tentials will be chan ging in time. Under these circumstances it 
is well known (e.g., ISachs&Wolfelll967l : iRees & Sciam3l 19681 : 



1 HEALPixURL site: http://healpix.jpl.nasa.gov 
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2 WMAP URL site: |http : // map . gsf c . nasa . gov| 
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iKofman & Starobinskiill 19851 : iMartinez-Gonzalez et al]| 19941) that 
CMB photons experience a gravitational blue/redshift, known as 
the Integrated Sachs- Wolfe (ISW) effect, i.e. 

Mr, n) 



w 2 r 
^ (n)= -J 



iir ■ 



dr 



(2) 



In this equation, ST denotes changes in the CMB blackbody bright- 
ness temperature, n refers to a given direction on the sky, and r de- 
notes comoving distance, which relates the time coordinate ; and 
the c onformal time n via t he rel ation dr\ = dt/a(t) = ±dr/c. 

ICrittenden & Turokl i ll 996b first realised that the presence of 
ISW could be separated from the CMB temperature anisotropics 
generated at the surface of last scattering by cross-correlating the 
measured CMB temperature map with density projections of the 
large scale structure distribution. The underlying idea is that if the 
initial metric perturbations generated during inflation are seeds of 
both, the potential and the density fluctuations, then these two fields 
should be correlated. This correlation can more clearly seen if both 
the ISW temperature anisotropies and the galaxy angular number 
density fluctuations are expressed in terms of the density contrast 
Fourier modes 6^. To perform these analysis both fields must be 
first decomposed on a spherical harmonic basis, 

ST> sw (h)IT = £ a' e s n y o „m, 8n g (h)/n g = £ a s (m .Y t , m (n), 

Ijn e,m 

(3) 

The multipole coefficients a' r s ™ , a s ,„, can then be written as (see, 
e.g.. ICooray||2002l; lHemandez-Monteagudoll2008l) 



«£T = H 



/rIV 



I 



dr j((kr) 



-3Cl,„Hl d(D(a)la), 



H)'(4. 



k 2 



f dk 

J (2jt. 



dr 



(4) 



J dr j ( (kr) W g (r) r 2 n g (r)b(r, k) D(r) 6 k j h g . 



(5) 

For each multipole I, m may take (2i + 1) different values, ranging 
from —I to t. n g (r) denotes the average comoving galaxy number 
density at the epoch corresponding to r; whereas h g denotes the 
average angular galaxy number density. The instrumental window 
function for galaxies at epochs given by r is provided by the func- 
tion W g (r). The functions jc(x) correspond to spherical Bessel func- 
tions of order t. In Equation 3, we have assumed that the galaxy 
field is a biased tracer of the matter density field, and thus the bias 
factor b(r,k) may a priori be a function of distance/time and/or 
scale k. If the galaxy and ISW maps are correlated, then the en- 
semble average of the product of the same (I, in) multipoles should 
yield a positive quantity, i.e., Cf isw = (a g f L(fi£? ')*) > 0, 



-gJSW 



R J k*dk P m (k) J dr, 

I 



Mkn) 



-3n„,Ho d(Dja) 



k 1 



dn 



dr 2 je(kr 2 ) r\ W g (r 2 ) n g (r 2 )b(r 2 ,k) D{r 2 ), 



(6) 



where P m {k) denotes the linear matter power spectrum at present. 
The auto power spectra are defined in a similar way, Cf = 
(afjafj'), with X, Y =ISW, CMB, g, etc, and the * symbol de- 
notes complex conjugate. The detection of this correlation between 



both maps (and hence the evidence for ISW) is, however, limited 
by the fact that the CMB is dominated by the signal generated at the 
Last Scattering Surface at z — 1, 050. In addition, the galaxy survey 
may be affected by shot noise, be incomplete or not probe the red- 
shift range where the ISW is generated. All these aspects must be 
taken into account when computing the theoretical signal-to-noise 
ratio (S/N) of the cross-correlation analysis. Since the signals de- 
scribed so far are written in multipole space, the S/N computation 
should be done naturally in multipole space as well. For each mul- 
tipole £, one readily finds that 



(") 2 



(Cf sw f 



x no.modes : 



( C f w fx(2{ + l)f sky 
Cf MB (C* + l/h g ) + (Cf' sw f 



(V) 

In this expression, the numerator contains the ISW - galaxy cross- 
correlation for multipole C, while the denominator refers to the vari- 
ance associated to its measurement. We next describe this equa- 
tion in detail. Provided that our theoretical description preserves 
isotropy, there is no m dependence in the power spectrum multi- 
poles. Furthermore, since there only exists one single universe to 
examine, the averages for each multipole ( are performed by con- 
sidering all corresponding (2t + 1) m modes (i.e., (...) = 1/(2^ + 
1) 2m •••)• However, if only a fraction of the sky f s fy is the subject of 
analysis, then the number of effective modes for a given multipole 
i shrinks by the same factor, i.e. (2£ + Y)f s k y . This explains why 
the variance for one single estimate of C s / S is divided in Eq.lfT} 
by (21 + 1)/^. Although this is a somewhat simplistic description 
of the impact of the sky mask on the S/N ratio, this approach has 
proved to capture the basics of the S/N degradation under incom- 
plete skies (see, e.g., ICabre et al.l2007l) . 

The variance associated to a single estimate of Cf ISW is given 
by the product of the CMB and galaxy auto power spectra plus 
the squared cross angular spectrum. The CMB angular power spec- 
trum is the addition of two (uncorrelated) components produced 
by the Surface of Last Scattering (SLS) and the ISW, i.e., Cf MB = 
Cf LS + C' ( sw . At the same time, provided that galaxies are discrete 
objects, the galaxy auto power spectrum contains a term account- 
ing for shot noise. This term should approximately obey Poissonian 
statistics and should be proportional to the average galaxy number 
density. However, since we are investigating relative fluctuations of 
the galaxy angular density field (i.e., 6n g (n)/h g and not Sn g (n)), the 
shot noise term in the power spectrum becomes inversely propor- 
tional to h g . 

From Eqs. d5|7t it is clear that the knowledge of the galaxy 
redshift distributions and the galaxy bias is required in order to 
provide an estimate for the ISW S/N ratio. If the bias is constant 
in space and time, then it is easy to show that it cancels out in 
Eq-0 provided the shot noise term is negligible. From the LOWZ 
and CM ASS LG samples we know the underlying spectroscopic 
redshift distribution, which has been obtained from a spectroscopic 
subsample, and is shown in Fig.lQJ). Our LG sample is practically 
insensitive to the onset of Dark Energy at z > 1, but it is able to 
trace its impact on the growth of potentials at z < 1, and this red- 
shift range becomes critical wh en evaluating the ISW contribution. 
iHemandez-M onteagudol J2008h found that the maximum contribu- 
tion to the ISW - galaxy cross-correlation was produced from red- 
shifts in the range z 6 [0.2, 1.0] when computing the ISW signal in 
redshift shells of different widths. The comoving number density of 
the LOWZ sample increases rapidly at low redshift, as displayed in 
Fig-GJi)- We consider only LGs placed at redshifts above a min- 
imum photometric redshift Zmm = 0.15, since our computations 
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Redshift z Multipole I Multipole I 



Figure 1. (a) Comoving LOWZ (red dotted line) and CMASS (green dashed line) LG number density as a function of redshift. (b) Angular power spectrum 
estimation of our LOWZ+CMASS sample. The blue dashed line corresponds to the shot noise Poisson term. The solid black line is the linear theory prediction 
for the Cf-s, adopting our fiducial ACDM cosmology model, and assuming the redshift distribution of our LOWZ+CMASS sample (solid line in panel (a)). 
The red line corresponds to the best fit to the observed angular power spectra (filled blue circles), which yields a constant bias estimate of b = 1.9. Error bars 
are approximated, for each multipole bin, as a fraction ^J2/(2f + 1)// S ky/A^ of the amplitude dictated by the red line, with M being the width of the multipole 
bin. See Appendix [B] for a description of an accurate error bar estimation, (c) Cumulative S/N ratio of the ISW - LG cross power spectrum below a given 
multipole £. Dotted, dashed and solid lines correspond to LOWZ, CMASS and LOWZ+CMASS samples, respectively. The latter converges to ~ 1.7 at I ^ 60. 



within the reference ACDM model show that, for a source popu- 
lation following the redshift distribution found for our LG samples, 
the total ISW S/N ratio possesses a local maximum around z m j„. We 
are effectively applying a cut in photometric redshift whose depar- 
ture from the spectroscopic one has negligible impact on the S/N 
estimates. Also, below this redshift, the LOWZ sample is known 
to become contamin ated by galaxies from the SDSS main sample, 
dPareikoetai]2013l) 

The cumulative ISW S/N ratio is greatly limited by our sky 
coverage: we must discard low Galactic latitude regions where 
WMAP CMB data may be contaminated, as well as those regions 
discarded by the SDSS DR8 survey footprint. Furthermore, since 
the cross-correlation analyses involving WMAP data require rela- 
tively low resolution maps with Nside=64, it becomes necessary to 
remove pixels which lie close to the borders as they drastically af- 
fect the estimated galaxy density. When combining both WMAP9's 
Kq85y9 mask and SDSS LOWZ+CMASS's effective mask, and af- 
ter dropping all mask pixels with mask value below 0.90, we have 
an effective sky fraction under analysis of / sky = 0.19. Finally, the 
linear constant bias for our LG sample is computed via a/ fit from 
the ratio of the measured LG angular power spectrum and the lin- 
ear theory prediction for the dark matter angular power spectrum, 
yielding a value of b = 1.88 ± 0.11, see Fig[]2b). 

We provide, in Appendix [B] a complete characterization of 
the angular power spectrum of the different galaxy samples, and 
analyze systematic residuals, study the consistency with the lin- 
ear Gaussian prediction and obtain constrains on local-type of non- 
Gaussianity. All these inputs allow the computation of the ampli- 
tude of the ISW - LG cross-correlation that we expect in the con- 
cordance ACDM model, and from it the differential (S/N) 2 per mul- 
tipole t as given by Eq.Q. In Fig.lQ]:) we show the integrated S/N 
ratio below a given multipole £, 



N 



(<£) = 



(8) 



The S/N converges at ~ 1.7 for t 
scale nature of the ISW signal. 



5 STATISTICAL METHODS 



60, as expected from the large 



In this section we outline the three different statistical methods we 
use in our cross-correlation analyses. One is based in real space, 
another in multipole space, and the third one in wavelet space. We 
end this section by applying these three analyses on mock ISW 
- LG data produced in the concordance ACDM model. HEALPi)0 
( Gors kTet al.l2005l) tools were used when implementing these anal- 
yses. In particular, unless otherwise indicated, all maps were down- 
graded to N siie = 64 resolution {9 pixel ~ 1°). 



5.1 The Angular Cross Power Spectrum (ACPS) 

As mentioned above, this method computes the average product of 
multipole coefficients from each (ISW and galaxy) map, i.e., 



f^g.CMB 



2C 



(9) 



In the presence of a sky mask, which excludes un-observed or 
contaminated regions, the estimates of the multipoles 2? , af™ 
become both biased and correlatecjfl as do the estimates of the 



C r ' -s. The im pact of the mask can either be corrected for (see 
iHivon et"al]|2002l) or included when comparing to theoretical ex- 
pectations. The computations of the covariance matrix for the 



4 HEALPixURL site: http://healpix.jpl.nasa.gov 

5 In the standard homogeneous and isotropic scenario (8^8 q) = 
(27r) 3 P m (k)6 D (k - q), and as a consequence for f s ty = 1 we have 
<of (al ,)*) = Cf Y S?„8 K ,; where S K . is the Kronecker delta (& K . = 1 if 

x f,m v t'jn ' £ 1,1' m,m" l,j v LJ 

i = j and Sfj = otherwise), i.e., different multipoles are uncorrected. 
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ACPS becomes necessary when estimating the statistical signifi- 
cance of the measured Cf CMB -s. In order to take into account in a 
realistic way the effects associated to the effective sky mask on the 
data, we follow a Monte Carlo (MC) approach. Ideally, one per- 
form MC mock catalogs for both CMB and LG components. In a 
first stage, however, we shall restrict our MC mocks to the CMB 
component. These MC CMB maps are uncorrelated to our LG cat- 
alog, and hence the variance of our ACPS estimates is given by 



qLO boss qCMB 
no.modes 



(10) 



This variance estimate approaches the denominator of Eq.((7j if 
Cf MB = Cf ,sw , (Cf w ? « Cf MB C s and the values of Cf are 
close to Cf G B0SS . In our ACDM cosmological model all these con- 
ditions should be fulfilled; if this is the case, running MC on the LG 
component should not appreciably change the variance estimate. 
Our approach can also be understood as a null hypothesis test, i.e., 
no existing correlation between our SDSS LG catalog and WMAP9 
data. In all cases, we ignore the / = 0, 1 multipoles (monopole and 
dipole). 

We may or may not average the estimates of Cf T into multi- 
pole bins centered around a given multipole lb, but in either case 
we attempt to add all the S/N within the relevant multipole range, 
which, according to Fig.lQJ), is chosen to lie in / 6 [2, Z„ rav ] with 



Int 



60. From each estimate C^' T , a corresponding rms value is 



computed from the MC simulations (cr« g ,r); using this rms value 

we calculate the following statistic that adds the S/N contribution 
from all multipoles: 



z 



C gJ 



S3 a q; 



(11) 



From our MC simulations we obtain an accurate estimate for the 
rms of this statistics a B[ as well, under the null hypothesis (i.e. 
(fie) = 0). The quoted total S/N is built from the estimate of Pi mas as 
the ratio 



S_ 

N 



(12) 



5.2 The Angular Correlation Function (ACF) 

This cross-correlation estimator, defined in real space, is computed 
from the subset of unmasked pixels from 



W(0) ■ 



Si 



je[B,g+de] 



2„ 



1 



(13) 



i,je[e,e+de] 

where the sum runs over the subset of pixels i, j that lie at a distance 
9. This method is related to the previous one: it can easily be shown 
that the full sky theoretical expectation for the ACF can be written 
in terms of the ACPS as 



w{0) : 



2t+ 1 



Z^ c r>Kcos,), 



(14) 



with P((x) the Legendre polynomial of order £. From this expres- 
sion, it is easy to see that ACF estimates at different angular separa- 
tions (9-s) will be highly correlated, since they are merely the sum 
of the same underlying Cf JSW -s with different weights (provided 
by the Legendre polynomials). The highest S/N is expected at zero 
lag (9 = 0), for which all multipoles contribute (Pf(V) = 1 for all 
€). As for the ACPS, in order to estimate the rms of the ACF esti- 
mates we first run MC on the CMB component as a null hypothesis 



R [degrees] 



1 .000 



0.100 r 



> 



7 0.010 




0.00 



Multipole I 

Figure 2. Multipole window function for Spherical Mexican Hat Wavelet 
of different sizes (R=20°, 10°, 5°and 1°) 



test. If an evidence for a cross-correlation is found, then MC on the 
galaxy template will be implemented. Both CMB and galaxy maps 
are normalized to have zero mean in the area under analysis. 

5.3 The Spherical Mexican Hat Wavelet (SMHW) 

The SMHWs form a family of filters defined, in real space, as 



wv(R, ft) 



1 



where R provides the filter scale, 
N(R) = R 



-(91 HI 



exp \-y 2 /(2R 2 )\ 



R 2 R 4 
1 + — + — , 
2 4 



(15) 



(16) 



and y is related to the angular distance 9 by y = 2 tan | . Each scale R 
defines a limited range of multipoles to which the filter is sensitive, 
as shown in Fig. {2}. Therefore, by first convolving the temperature 
and galaxy maps with these filters and then computing the cross- 
correlation of the resulting maps, one obtains constraints on cross- 
correlation at different angular scales. One should, however, keep 
in mind that these angular ranges are not disjoint, and hence results 
for different values of R will a priori be cor related. For the firs t 
implementation of SMHW in ISW studies, see lVielva et alj J2006h . 
while further insight in the SMH W properties can be found in e.g., 
iMartfnez-Gonzalez et alj d2002l) . This implementation is relatively 
simple: multipole coefficients for the wavelet at each scale R are 
first computed (wv (m {R)). Given the azimuthally symmetric form 
of Eg.lll5b. only m = multipoles contribute and the convolution 
is performed by simply multiplying the multipole estimates of , 
tf m B b v me corresponding wavelet counterpart wv fj „ =0 (R). Maps 
are inverted back onto real space, and the zero lag cross-correlation 
is computed as the average product outside the mask, by 



n(R) = 



Z,ST(R,n,)S g (R,ii,) 



(17) 



where s(R, ft) denotes the convolved version of map s(h) with 
SMHW at scale R, and the sum runs only through pixels outside 
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the maslfl As in previous cases, MC CMB mock maps were pro- 
duced and processed exactly as WMAP9 data, and the outcome of 
these tests provides the null case probability distribution function 
for the IT(i?)-s, to which observed data results are to be compared. 
In particular, those simulations provide the rms for each R scale 

As for the ACPS, we build the statistic 



z 



(18) 



where the indexes i and j run for the different filter scales R con- 
sidered. The different terms in the sum of the right hand side of this 
equation are correlated, although the MC simulations account for 
those correlations and are able to provide an estimate of the rms for 
this statistics erg . After adding the contribution from all scales, we 
quote the total S/N ratio obtained with this statistics, i.e., 



S_ 

N 



(19) 



I'm 



where the subindex "nf " denotes the total number of filters consid- 
ered. 



6 RESULTS 

In this section we present the results obtained with the three dif- 
ferent cross-correlation methods described above, and is divided 
in two sub-sections. In the first subsection, we build mock galaxy 
catalogs following the SDSS-III/BOSS LG clustering and redshift 
distribution properties (according to our reference ACDM cosmol- 
ogy). These LG mock catalogs have shot noise levels identical to 
those in real data. Each of those galaxy mocks is correlated to a sim- 
ulated Gaussian ISW map following the angular power spectrum 
computed in the same cosmological model, and each simulated 
ISW map is added to a Gaussian CMB component according to 
the angular power spectrum of the CMB anisotropies generated at 
the surface of Last Scattering. Thus, we have two sets of simulated 
maps, one for galaxies and one for CMB. Real sky masks are ap- 
plied when feeding these simulated maps in our analysis pipeline. 
The second subsection presents the results obtained from observed 
data, under different levels of corrections for systematics, and are to 
be compared to the results of the previous subsection. The impact 
of other artifacts associated to the effect of the mask under use is 
also addressed in this section, while the description of the attempts 
to correct for potential systematics is provided in Appendix [A] As 
mentioned above, most of the cross-correlation analyses were per- 
formed at the more easily manageable map resolution of Nside=64. 
However, some consistency tests were run under Nside = 128. 



6.1 Results from ACDM simulations 

Fig. [3] displays what should be expected for our cross-correlation 
analysis in a typical concordance ACDM simulation. One random 



6 In practice, the sum is taken only through those pixels for which the con- 
volved mask (with the corresponding wavelet) is above a given threshold. 
We show results for a threshold set equal to zero, but the results change 
only by a few percent under a threshold in the range [0, 0.5]. Our threshold 
choice is motivated by simulations, which show that higher S/N is recovered 
when including all partially masked pixels in the analysis. 



LG mock catalog was analyzed in conjunction with a correspond- 
ingly simulated (and correlated) CMB map, under the same con- 
ditions of sky coverage as the observed data. The left panel shows 
the results for the ACPS analysis: red circles denote positive values 
for the ACPS multipoles, while blue squares correspond to negative 
values. This convention is followed in the other panels. The dashed 
and solid lines correspond to the 1 and 2 cr levels, respectively, ob- 
tained after running 10,000 MC CMB simulations under the null 
hypothesis, i.e., our LG mock catalog is kept fixed and CMB skies 
are generated and cross-correlated to it. This is also the procedure 
used when interpreting the results from the SMHW and ACF tests. 
Ax 2 computation on our ACPS estimates w.r.t the null hypothesis 
(for which all Cf isw = 0) yields^ 2 = 28.1. This is obtained after 
grouping all multipoles in the 32 different bins shown in the fig- 
ure: about 66 % of the null simulations yield higher values of x 2 , 
so this test does not seem to be too sensitive to the ISW-induced 
correlation. The /3 statistic defined in Eg. ill It yields a S/N of 2.19 
(according to Eq. 1 1 2b . which reflects the presence of at least five 
positive points at the 2 cr level in the left panel of Fig. [3] together 
with a majority of positive values of the cross-spectrum estimates. 

In the middle panel of Fig. [3] we display the results for the 
SMHW coefficients. On all scales the coefficients are positive, and, 
on the largest scales, they cross the 2 cr level. In this case, for the 
nine filter scales under consideration we obtain^- 2 = 13.96 (about 
12 % of the simulations provided higher values of x 2 )- The S/N ob- 
tained from the /j nf statistics, introduced in Eq.dl9|l, yields 1 .65, and 
4.82 % of the 10,000 MC null simulations provided higher values 
of y6„f. Finally, the ACF analysis is based upon only 100 MC null 
simulations, which are nevertheless sufficient to define the 1 and 
2 cr levels. In this case, the ACF goes above the 2 cr level in the low 
9 range (S/N=2.26 at 8 = 0°, see right panel of Fig.[3}- 

The values obtained for the S/N in the ACPS and SMHW anal- 
ysis are in good agreement, within 0.5 cr, with the theoretical pre- 
diction of S/N^ 1.67 quoted in Section|4] After running 10,000 ad- 
ditional MC simulations containing correlated CMB maps and LG 
mocks, we obtain average S/N values of 1 .48 and 1 .44 for the ACPS 
and SMHW algorithms, respectively. Likewise, the average zero 
lag S/N of 100 mock estimates of the ACF amounts to 1.48, which 
is in excellent agreement with the other two algorithms. These S/N 
average values are within 14 % of the (approximate) theoretical es- 
timate of 1.67 obtained from Eq.[8] 

One further issue to address is to what extent the corrections 
applied on the LG templates would affect the S/N of our cross- 
correlation estimators. In Appendix |A] we describe the procedure 
by which the LG number density is corrected from possible correla- 
tions to known potential systematics (such as stars, sky emission or 
dust extinction). We use the difference between the raw LG density 
map and the one accounting for all possible corrections (namely 
star density, mask pixel value, seeing, sky emission, dust extinction 
and airmass; see Appendix[A}. This difference map (Fig. lA5t con- 
stitutes our template for the systematic impact on our mock galaxy 
templates. We run again the same set of MC simulations accounting 
for the correlation between the galaxy mocks and the CMB reali- 
sations, but after adding our template for systematics to each sim- 
ulated galaxy number density map. In this new set of simulations, 
the average values for the S/N statistics amount to 1.20 and 1.19 
for the ACPS and SMWH approaches, respectively; while the av- 
erage zero lag amplitude of 100 MC mocks of the ACF yield a S/N 
of 1.03. This exercise demonstrates that the average impact of the 
systematics that our analysis has identified and corrected amounts 
to a 17 %-30 % decrease of the average S/N. Although significant, 
these changes are not as critical as those induced in the auto power 
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Figure 3. Results from a single ACDM simulation; in all panels red (blue) symbols denote positive (negative) values. Left panel: ACPS results for / e [2, 60]. 
Solid and dashed lines display the 1, 2tr confidence level after 10,000 MC CMB simulations under the null hypothesis. When combining all multipoles, we 
find an integrated significance at the S/N 2.19. Middle panel: Results for the SMHW test after considering scales of R = 3.3°, 4.2°, 5.0°, 6.6°, 8.3°, 10.0°, 
12.5°, 15.0° and 17.5°. Again, solid and dashed lines display the 1, 2-cr confidence level after 10,000 MC CMB simulations under the null hypothesis. When 
gathering information at all scales, we obtain S/N ^ 1.65. Right panel: Results of the ACF analysis. Dashed and solid lines have same meaning as in previous 
panels, although these are built upon only 100 MC simulations. At the low 6 range the ACF lies around the 2.2 tr significance level, in good agreement with 
the other two statistical methods. 




Figure 4. Results from cross-correlating WMAP-9year V-band with LG templates constructed from LOWZ and CMASS samples after correcting for stars, 
mask value, seeing, sky emission, dust extinction and airmass. In all panels red (blue) symbols denote positive (negative) values, and dashed and solid lines 
display the 1 <x and 2 tr confidence levels, respectively. Panel organization is identical to Fig. [3] In this case, the S/N ratio obtain from the ACPS, SMHW and 
ACF methods lie at the 1 .62 - 1 .67 tr level, in good agreement with the ACDM predictions. 



spectra of the galaxy samples at the low I regime, (see Figs. IA3I 
IA4l in AppendixlAt. 



6.2 Results from observations 

In agreement with the tests described above, an investigation of the 
impact of the correction for systematics in the recovery of the S/N 
of the LG - ISW cross-correlation shows that, on the observed data 
set, the corrections for potential systematics described in Appendix 
|A]have an effect at the 23 % - 30 % level of the total S/N. 

In Fig. [4] we display the results, under Nside=64, for the cor- 
rected galaxy template. In this case, the,^ 2 test yields 38.1 for 32 
d.o.f., i.e., not able to rule out by itself the null hypothesis (more 



than 21 % of the null-case MC simulations provided higher,^ 2 val- 
ues). In terms of the /? statistic for the ACPS, the S/N becomes 
1.62 (only 4.85 % of the mock simulations provided larger values 
of this statistic). For the SMHWs,^ 2 = 9.34; this test is again less 
conclusive than the fi n{ statistics, which yields S/N= 1.67. Finally, 
the zero-lag amplitude of the ACF analysis results in S/N= 1.69, 
close to the other two estimates. Had no correction for systematics 
been applied, the attained S/N levels would have been considerably 
lower. For instance, for the ACPS and SMHW p statistics, the S/N 
values fall to 1.24 and 1.16, respectively, while the zero-lag ampli- 
tude of the ACF yields S/N= 1.27. A comparison between Fig. [4] 
and Fig. [5] shows that the artifacts present in the raw maps only 
slightly bias the cross-correlation estimates, while widening con- 
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Multipole t Wavelet scale [degrees] d [degrees] 

Figure 5. Same as Fig. [4] but after using the raw LG templates, before applying any correction. The overall S/N level decreases significantly (~ 0.5 cr) when 
compared to the corrected one shown in Fig. @] 



SMHW" product map at K=6.6 deg 
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Figure 6. Mollweide projection in Galactic coordinates of the wavelet prod- 
uct map at a scale of R= 6.6° from WMAP-9yr V and LOWZ+CMASS 
data. Units are in the rms of this product map, and the central longitude 
corresponds to I = 0°. 

siderably the allowed area for the 1 and 2 cr confidence levels (i.e., 
increasing the effective error bars of the cross-correlation analysis). 

The wavelet product maps provide a visual description of the 
areas where the cross-correlation between the galaxy and CMB 
maps is built. In Fig. {3} we display the product of the wavelet 
convolved maps at the scale of 6.6°, in units of the rms of the same 
map. The excess of positive spots over negative ones gives rise to 
the ~ 2 cr excursion at that scale in the middle panel of Fig. (|4j. 

We next explore how the recovered S/N levels depend on the 
mask choice. In our standard mask, we consider for analysis pix- 
els with mask values at Nside=64 ranging from 0.9 up to 1. We 
have checked that making the mask equal to one in all those pix- 
els introduces negligible changes (below the percent level in S/N). 
Since downgrading the original mask to Nside=64 introduces bor- 
der effects that force a ~ 16 % decrease of the effective area under 
analysis, we conduct a single ACPS and SMHW test at Nside=128 
resolution. At this resolution, the allowed sky fraction where border 
effects are kept under control increases to / sky a 0.23. If S/N truly 
scales as oc ^/ sky , then one would expect a ~ 7-10 % improvement 
of the S/N. In practice, for the corrected LG template we obtain, 



after 1,000 MC null simulations, S/N values of 1.97 and 1.77 for 
the ACPS and SMHW tests, respectively. These values constitute 
an increase of 0.35 cr and 0. 10 cr w.r.t the Nside=64 case. While 
the SMHW S/N increase is in perfect agreement with the predicted 
/sky scaling, the ACPS one is considerably higher (~ 0.2 it) than 
expected, although in the appropriate direction. 

Even under Nside=128, our results cannot rule out by them- 
selves an Einstein-de Sitter scenario at the 2 cr level. However, the 
recovered cross-correlation level is in excellent agreement with the- 
oretical expectations based upon the concordance ACDM cosmol- 
ogy- 



7 DISCUSSION AND CONCLUSIONS 

After implementing three different statistical approaches to mea- 
sure the cross-correlation between the CMB and our SDSS- 
III/BOSS luminous galaxy template, we have found good consis- 
tency among their outputs both on simulated and observed data. 
The mock simulations have shown levels of S/N close to the theo- 
retical prediction of S/N^ 1.67. In particular, an ensemble of coher- 
ent simulations have provided average levels of correlation within 
14% of that theoretical estimate, while the analysis of a single 
random mock realization also yields consistent significance levels 
within half cr from the expected value. The impact that inaccuracies 
in the systematic removal may have on the obtained S/N is found to 
be at the ~ 20-30 % level for our mock galaxy samples, and hence 
it becomes of relevance when comparing accurate model predic- 
tions obtained from numerical simulations with analysis obtained 
from observational data. 

The solid consistency among the three different statisti- 
cal methods provides evidence for the robustness to the cross- 
correlation results involving WMAP-9 yr data and SDSS-III/BOSS 
data. For all three ACPS, SMHW and ACF analyses, the S/N lev- 
els are found to lie between 1.62 - 1.69 cr. Furthermore, we have 
found that not correcting for the impact of systematics on the mea- 
sured LG number density does change the low £ angular auto power 
spectrum estimates; this translates into a lost of S/N at the ~ 0.5 cr 
level. 

Our results on the observational data are in good agreement 
with the ACDM predictions since they lie within a half cr from the- 
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oretical expectations, although they do not exclude an Einstein de 
Sitter scenario by themselves. Our results are in clear tension, how- 
ever, with results obtained from earlier SDSS data sets. Somewhat 
peculiarly, the significance of ISW detection in SDSS data sets has 
tended to decrease as the data set has gotten larger. iFosalba et alj 
d2003h reported a 3c ISW detecti on using an early release of 
the SDSS main galaxy sample and Cabr e et al] d2006h reported 
4.7 cr detection using SDSS DR5 data. More recent claims from 
Ho et alJd2008h (at - 2.46 cr utilizing high redshift SDSS LRGs) or 
Giannantonio et al. j2012l) (~ 2.5 cr adopting the LRG DR7 Sloan 
sample) are more modest, but still fall ~ 0.7 <r above the theoretical 
predictions inferred from our photoz SDSS DR8 sample, which is 
a deeper and wider tracer of the lar ge scale gravitational potentials. 
Actuallv. lGiannantonio et al] d2012h point out that their results with 
LRGs from SDSS DR7 are about 1.3 cr above the theoretical S/N 
ratio expected for that surv ey. On the other hand, th e 'SDSS' and 
'2slaq' samples analyzed in Sa wangwit et alj d2010h are similar in 
depth to our LOWZ and CM ASS samples and their ACF measure- 
ments for those two sampl es agree well with A CDM predictions. 
In a parallel analysis Giannanto nio et all d2013l) analyzed a DR8 
LRG sample similar to our own, and found a similar significance 
to our own (1.8) and this result is in good agreement with ACDM 
predictions. 

The reason for decrease in the ISW significance obtained with 
DR8 data is unclear. While in previous cases it was speculated that 
the excess anisotropy power on large scales may be responsible for 
an increased level of evidence for the ISW, in our case we meet 
the opposite situation, i.e., the presence of artifact systematics in- 
creases the level of anisotropy on large scales, and consequently 
also increases the errors/uncertainties in the recovery of the cross- 
correlation coefficients in that angular regime. The amplitude of the 
cross-correlation is not so critically modified, and thus we conclude 
that, at least for our SDSS DR8 data set, the presence of artifacts in 
the galaxy template tends to worsen the ISW S/N constraints. 

In this analysis, we use ACPS, SMHW and ACF measure- 
ments and find similar results with each technique. One possible 
drawback of the use of the ACF as unique cross-correlation anal- 
ysis is that it is sensitive to all scales present in the map, in par- 
ticular to the small ones which, while not carrying much ISW in- 
formatio n, are the most subject to point s ource contamination. For 
example, Hernandez-Monteagudc] d20ld) found that, when cross- 
correlating a radio survey to WMAP data, the ACF provided much 
higher levels of cross-correlation than the ACPS, and the signifi- 
cance level ( i) was correlated to the flux threshold applied to radio 
sources and (ii) restricted to the smallest scales. This example il- 
lustrates the utility of implementing simultaneously Fourier- and 
real-space based algorithms in ISW studies, although it does not 
imply that point sources must necessarily be biasing all previous 
analyses. In particular, in our galaxy sample, we find no evidence 
for significant point source emission in WMAP-9yr sky maps as 
suggested by the fact that our ACF significance is almost identical 
to the ACPS significance. 

Another possible reason for the mismatch with respect to some 
of the previous works may be associated to the particular statistic 
used to infer the ISW significance level and its dependence on the 
assumed model describing the clustering of the LGs. In our case 
we are quoting deviations wrt the null hypothesis and our quotes 
are insensitive to the theoretical model describing the LG angula r 
power spectrum. In other works (e.g.. iGiannantonio et al.ll2012h . 
a fit to a theoretical prediction for the power spectrum/correlation 
function of the measurements is performed, and different assump- 



tions in such theoretical modeling may give rise to different S/N 
quotes. 

In our study, we make use of all the relevant information avail- 
able in the SDSS DR8 to not only provide a clean tracer of the 
potential wells at intermediate redshifts, but also to make precise 
predictions on what should be seen with this template in an ISW - 
cross-correlation study in the ACDM cosmology. The use of Gaus- 
sian simulations combined with Poisson sampling should further 
refine the predictions for the level of ISW signature on the large 
scale structure. Although our results are in good agreement with 
expectations, they do not provide conclusive evidence for the pres- 
ence of ISW. 

In practical terms, we conclude that the study of the ISW ef- 
fect becomes limited by the accuracy to which a given galaxy sur- 
vey can characterize its source distribution. Uncertainties not only 
in the galaxy bias, redshift distribution and evolution, but also the 
impact of sources of systematic uncertainty on large-scales, are cur- 
rently higher than the uncertainties in the CMB maps, and hence 
they control the prospects for the ISW characterization. This is- 
sue reaches greater relevance when different surveys (each of them 
with their own uncertai nties) are comb i ned in a single ISW cross - 
correlation study, (e.g. iHo et al.ll200ll ; IGiannantonio et al.ll2012h . 
The future of ISW science does not seem to reside so much on the 
CMB data, but on the side of an accurate cartography of the large 
scale structure of the universe. 
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APPENDIX A: SYSTEMATICS IN THE LG DATA 

It is well known that the measured SDSS DR8 LG density on the 
sky is affec ted by several sys t ematics (see the study for the CMASS 
sample of iRoss et alJl20lH IHo et alj[2oT3) . The impact of stars 
is known to be the most relevant, although in this Section we at- 
tempt to correct for the biases introduced by other possible system- 
atics such as the effective value of the mask in each pixel, the see- 
ing, the sky emission, the dust extinction and the airmass. The star 
density template was built upon a star catalogue restricting to the 
imod magnitude range [17.5, 19.9]. The mask used in our analysis 
only considers pixels at iVside = 128 resolution having mask value 
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Figure Al. Scaling of the LG number density versus different potential systematics, namely star density, mask pixel amplitude and seeing (first to third 
column). Top row panels provide the fraction of the observed area versus the amplitude of the potential systematic, such that the integral below the histograms 
should equal unity. The middle and bottom rows refer to the CMASS and LOWZ samples, respectively. Black (blue) circles display the scaling before (after) 
correction (see text). Red solid line displays a spline fit to the measured scaling of the LG number density wrt each potential systematic. 



above 0.85. There exist records of the average seeing at each ob- 
served pixel (given in arcsecs), and of the average sky emission in 
the ;-band (given in nanomaggies arcsec -2 ), together with the dust- 
in duced extinction in t he r-band, A r , derived from the dust maps 
of lSchleeel et alj dl998h . Finally, we also considered the airmass as 
potential source of LG number density modulation. 

In our approach we investigate the scaling of the LG density 
w.r.t quantities that are possible sources of systematics, separately 
for the LOWZ and CMASS samples. When studying the signifi- 
cance of this scaling it is assumed that LGs are Poisson distributed, 
i.e., we neglect their intrinsic clustering. This approach allows us to 
assign error bars to each scaling in a cosmology-independent way, 
which are however too small, and this must be present when inter- 
preting the scalings. 

When building the raw CMASS galaxy density template, we 
weight each galaxy by its star-galaxy separation parameter (p sg ), 
which provides the pro bability of a given object to be mistaken by 
a star, (Ros s et alj201ll) . This parameter was, however, ignored for 
the LOWZ sample, since in this case there is significantly smaller 
chance for confusing a galaxy with a star. In Fig. IA1I we study 
the variation of LG density w.r.t star density, the mask pixel value 
and the seeing. The top row shows the relative fraction of the ob- 
served sky in each bin of each potential systematic (star density, 
mask value and seeing value), in such a way that the integral under 
the histogram equals unity. The middle (bottom) row corresponds 



to the CMASS (LOWZ) sample. In the middle and bottom panel 
rows, the black filled circles denote the scaling of LG density ver- 
sus the corresponding potential systematic. The red solid line, in 
each panel, displays the spline fit to that scaling, and the blue filled 
circles show the result after correcting the LG density field by the 
scaling provided by the spline fit. The LG density field corrected 



for the "X" potential systematic is given by rig c (n) in 



" g (n) 
S[X(n)]' 



(Al) 



where n g (n) denotes the raw LG number density. The vector n de- 
notes a given direction/pixel on the sky. The symbol S [X(n)] de- 
notes the spline fit of the LG number density - systematic "X" 
scaling evaluated at the value X(n). The normalization for this cor- 
rection is such that the average LG number density before and 
after correction are identical. As mentioned above, error bars as- 
sume Poissonian statistics. This figure format is identical to that of 
Fig. lA2l where the correction for sky emission, dust extinction and 
airmass are considered. 

Each LG density map is corrected for each systematic in suc- 
cession as displayed from left to right in Figs. dAl IA2t . The order 
at which each potential systematic is corrected should not matter 
if the effects are independent. However, stars and dust extinction 
are strongly correlated to Galactic latitude in a similar way, so this 
analysis requires more caution. From Figs. (1A1IA21 it is clear that 
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Figure A2. Same as in Fig. lAll but referred to sky emission, dust extinction and airmass. 



stars seem to be significantly modulating the LG number density, 
since for both LOWZ and CMASS samples the LG number density 
decreases rapidly with star density (although admittedly the effect 
is slightly stronger for the LOWZ sample). In either case, the cor- 
rected number density seems not to show any remarkable depen- 
dence on star density on more than 99 % of the observed area (see 
blue circles below ~ 900 deg~ 2 ). As shown in the middle panels, 
the LG number density appears to steadily decrease for increasing 
values of the mask, suggesting a too conservative estimate of the 
amount of area masked out around holes present in the footprint. 
This is common for LOWZ and CMASS samples. When studying 
this scaling for the degraded Nside=64, this scaling is flipped, and 
at low mask pixel values one finds lower galaxy densities, pointing 
to a clear footprint border effect. The seeing, conversely, does not 
appear to significantly modify the LG number density: only for the 
CMASS sample are there any hints of a trend by which less LGs 
are identified in regions with worse seeing, although this induces a 
small correction in the resulting LG number densities. 

The left column in Fig. lA2l shows that the sky emission modu- 
lates the CMASS LG number density in a non-intuitive way: the 
CMASS number density appears to increase with sky emission. 
However, this trend is based upon a relatively small fraction of the 
observed area (< 8 %), so, if real, it should have a relatively mod- 
est impact. The LOWZ sample does not show evidence for such a 
behavior. On the other hand, dust appears to influence the LOWZ 
number density but not the CMASS one (see middle column in 
Fig. lA2t . Since dust and stars are spatially correlated on the sky, we 
explore the scaling of LG density w.r.t dust extinction before and 



after correcting for the star density: we find no evidence for any 
significant change in the LG density versus dust extinction scaling. 
Only the LOWZ sample scaling suggests that dust may induce con- 
fusion in the identification of LOWZ LGs, since the relation of the 
CMASS number density with dust extinction is essentially flat. Fi- 
nally, the airmass seems not to introduce a significant trend in the 
LG number density, as shown by the right column of Fig. lA2l 

In Figs. J A3 1 IA4b we display the induced changes in the re- 
sulting angular power spectra. These spectra are computed after 
correcting for the bia s induced by the s ky mask by means of the 
MASTER approach Jffivon et alj|2002h . Since we shall perform 
cross-correlation studies with CMB maps, and in these maps the 
dipole has no cosmological information, we choose to remove the 
residual dipole in the effective area in both CMB and LG sky maps 
using the remove_dipole routine of the HEALPix distribution. In 
those plots, the correction for stars introduces the largest changes in 
the estimated angular power spectra, since the quadrupole (I = 2) 
and the multipole band centred at £ = 5 drop by factors of ~ 20 
and 6, respectively, for the CMASS sample (changes for the LowZ 
sample are much more modest). Practically all corrections leave all 
band power above I = 5 effectively unchanged, while the first two 
band power spectra display a more unstable behavior. After star 
correcting the CMASS sample, the band power centred at multi- 
poles I = 2,5 typically shift between 4 x 10~ 6 and 3 x 10~ 5 . For 
the LOWZ sample the band power centred at I = 5 barely shifts 
from ~ 1.2 x 10~ 5 , while the one centred on the quadrupole (t = 2) 
is negative until the dust extinction correction is applied (Negative 
cases are denoted by crosses rather than filled circles in those fig- 
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Figure A3. Impact of the correction for stars, mask pixel value and seeing on the estimated angular power spectra. Black and red circles correspond to angular 
power spectrum estimates before and after the correction, respectively. Crosses denote negative values, which may arise after correcting for the mask in the 
MASTER algorithm. 



ures). The MASTER algorithm may provide negative outputs of 
auto-power band spectra for cases where these are very poorly con- 
strained (as is the case here). Corrections are applied sequentially, 
so that the last correction is made upon a map previously corrected 
for all previous possible systematics. 

The total changes in our template for angular galaxy density 
can be seen in Fig. JA5b . Clearly visible is a low Galactic latitude 
increase in the galaxy density (a compensation for the star-induced 
galaxy obscuration). Also visible are stripes corresponding to the 
survey scanning strategy. 



APPENDIX B: THE ANGULAR POWER SPECTRA OF 
THE LG DENSITY TEMPLATES 

In this Section we address the computation of the angular power 
spectrum of the LOWZ and CMASS LG samples after conducting 
the corrections for systematics. We use the MASTER algorithm to 
correct for the bias induced by the joint sky mask (respecting only 
a ~ 23 % of the sky) on the measured angular power spectrum. The 
outputs of the MASTER algorithm in each of our galaxy samples 
are displayed by red crosses in the three panels of Fig. (1A61 for 
LOWZ, CMASS and LOWZ+CMASS samples, from left to right, 
respectively. The entire multipole range is divided in different bins, 
in which MASTER provides an estimate of the band power spec- 
trum. 



Using the redshift distribution obtained for each of the galaxy 
samples displayed in the left panel of Fig. Q}, we compute the lin- 
ear prediction for the angular power spectrum of the angular den- 
sity contrast. We next fit a constant bias for each galaxy sample: in 
order to minimize the impact of non-linear effects, we restrict our- 
selves to the band power spectra contained in the multipole range 
I 6 [2, 100]. Each band power estimate is initially assigned an er- 
ror equal to C t ^2/(21 + l)// sky /A£, where C c is the MASTER out- 
put for the band power estimate, C is the bin central multipole and 
hi its width. With these error estimates, we perform a x 2 fit to a 
constant bias relating the observed power spectra and linear the- 
ory predictions in our concordance WMAP-9yr ACDM model. For 
the LOWZ, CMASS and LOWZ+CMASS samples, we obtain con- 
stant bias values of b = 1.98 ± 0.1 1, 2.08 ± 0.14 and 1.88 ±0.11, 
respectively. As predi cted by the scalin g of the bias versus comov- 
ing number density of Nuz a et alj | |2012|) . the measured bias values 
decrease for the sample having higher number density. Although 
such scaling is referred to the bias computed in space, in our case 
we find that the joint LOWZ+CMASS sample, having higher an- 
gular density than the other two, yields a slightly smaller bias value 
as computed from the ratio of angular power spectra. It is neverthe- 
less worth noting that the bias estimates do not deviate from each 
other more than 2 a. 

With these values for the LG bias we next perform 1,000 
MC full-sky Gaussian simulations of galaxy density, which are re- 
normalized to have the observed average galaxy number density of 
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Figure A4. Same as in Fig. lA3l but referred to sky emission, dust extinction and airmass. 



each sample. These full-sky Gaussian simulations follow an input 
angular power spectrum and are produced with the HEALPix tools, 
which introduce some smearing due to the pixel window function in 
the simulated maps. We further run Poissonian realizations of the 
expected number of galaxies in each map pixel, hence providing 
our galaxy density templates the corresponding level of shot noise 
predicted by Poissonian statistics. Each full sky map is multiplied 
by the effective sky mask used in observed data analysis, and the 
monopole and dipole are removed from the remaining sky fraction. 

The resulting map is then analyzed with the MASTER algo- 
rithm, which produces estimates of the angular power spectrum 
from our surviving sky fraction. After correcting for the pixel win- 
dow function and the shot noise term (which equals 1 /», with n the 
average galaxy number density in each sample), we obtain the av- 
erage recovered angular band power spectra estimates displayed by 
filled black circles in Fig. dA6t . These symbols closely follow the 
linear prediction of the angular power spectra for each sample, after 
being boosted by the square of the corresponding bias factor (see 
solid black lines in each panel). Only the first band power estimate 
(centred at C = 2) clearly deviates from the linear prediction: this is 
a consequence of removing the dipole within our effective sky cov- 
erage. Since f& y < 1, the dipole is coupled to immediately higher 
multipoles, and removing the dipole in our effective sky biases low 
the first band power estimate. The rms of each band power estimate 
around the average value is displayed as the effective error bar as- 
signed to the real band power estimates, depicted by red crosses in 
Fig. dA6t . Finally, shot noise levels are depicted, for each sample, 
by horizontal dashed lines. 



We next address the question on how consistent is the data 
with the linearly biased ACDM linear theory prediction for each of 
the three galaxy surveys. For that, we use the 1,000 MC simulations 
described above to construct a covariance matrix between the dif- 
ferent band power spectra estimates. This allows conducting a x 2 
test for the measured angular band power spectra for each galaxy 
sample on the (biased) linear theory predictions. When restricting 
to multipoles below £ = 100, non-linear effects are not important 
and we obtain, for the 21 band power spectra estimates, x 2 val- 
ues of 13, 30 and 26 for the LOWZ, CMASS and LOWZ+CMASS 
samples, respectively. For these galaxy samples, we find that 91 %, 
12% and 26% of the MC simulations yield higher^ 2 values than 
those produced by the observed data, i.e., the measured angular 
power spectra are compatible with the linearly biased model predic- 
tions. However, when a wider multipole range is considered, then 
non-linear effects appear and the x 2 test yields worse results: for 
all band power spectra below t = 154, 14 %, 0.3 % and 29 % of the 
MC simulations provide higher^ 2 values. Furthermore, when look- 
ing at the entire range C e [2, 205], those ratios become 0%, 0% 
and 3.5 % for the LOWZ, CMASS and LOWZ+CMASS samples, 
respectively. Clearly, by including smaller scales one considerably 
worsens the fit. Out of the three samples, the most noisy seems to be 
the CMASS sample, with typically larger^ 2 values, while lower x 2 
values are produced by the LOWZ sample when remaining below 
£ = 100. 

The observed data points in Fig. jA6) are overall in good 
agreement with the corresponding linear theory predictions in our 
Gaussian concordance ACDM model, pointing to a low level of 
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Figure A5. Mollweide projection in Galactic coordinates of the effective correction applied to the combined LOWZ - CMASS galaxy templates, in units of 
deg -2 , with central Galactic longitude at I = 0. 
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Figure A6. Angular power spectrum estimates for the LOWZ (left panel), CMASS (central panel) and LOWZ+CMASS (right panel) samples. Observed data 
are displayed by red squares, while the average of the output of 1,000 MC simulations are provided by filled black circles. The best-fit linearly biased ACDM 
prediction is provided by the black solid line; it constitutes the input for the MC simulations. Finally, the shot noise contribution is indicated by horizontal 
solid lines in each case. 
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primordial non-Gaussianity. We next restrict our analysis to the 
combined LOWZ+CMASS s ample. We i mplem ent the local non- 
Gaussian b ias model u s ed inlRoss et alj J2013[). which w as itself 
inspired by iDalal et all d2008h ; iMatarrese & Verdel d2008l) . In this 
model, the effective bias is the addition of a (constant) Gaussian 
bias to a non-Gaussian contribution, 



b+{b-\)f™ 



(Bl) 



k 2 T(k)c 2 ' 

where b constitutes the Gaussian bias, £1 is the total matter critical 
density parameter, Hq is the Hubble constant and S c (z) oc 1 /D(z) is 
the spherical collapse overdensity, with D(z) the linear matter den- 
sity growth factor. The function T(k) is approximated by the dark 
matter transfer function, k is the wavenumber in units of Mpc~ 1 and 
f^f is the local non-Gaussianity parameter. We compute estimates 
of the angular power spectrum for the LOWZ+CMASS sample for 
different values of the Gaussian bias (b e [1.70,2.30] and the lo- 
cal non-Gaussianity parameter C/j!j£ al 6 [-300,300]). We perform 
a direct computation of the posterior probability for the bias and 
/jljj™ 1 after assuming Gaussian priors for the bias and assigning them 
rms values equal to those found in the^ 2 test described above. We 
find, however, that our constraints are practically independent of 
the choice of these priors, and since uncertainties in the amplitude 
of the anisotropy normalization can also be attached to this parame- 
ter, our results are insensitive to this normalization. We restrict our 
analysis to the multipole range C e [4, 100], and ignore the con- 
tribution of systematics to the covariance, since even for the band 
power centred on £ = 5 the extra error induced by uncertainties 
in the correction procedure is in all but one cases negligible when 
compared to cosmic variance. For the covariance matrix of the Q-s 
we compute the covariance of the 1,000 MC Gaussian simulations 
of our ACMD reference model, and scale it to the amplitude of the 
theoretical angular power spectrum of each model under consider- 
ation, 



at 95 % C.L.), although in that case the method of handling and 
marginalizing for syste matics is different. Although the works of 
IXia et alJ 120103 . 1201 lh may have not considered explicitly the sys- 
tematics under study in our work, they have excluded very low mul- 
tipoles from their analysis, and yet they find som e weak evidenc e 
for positive /^ cal , (/^ cal = 48 ± 20 at 68%C.L..IXia et al 11201 lh . 
Using a similar set of galaxies, Gian nantonio et alj ( 120131) found 
results consistent with ours, (-90 < /^ cal < 120 at 95 %C.L. for 
their LRG only sample). 

We finally address the issue of the North to South asymme- 
try. Our mock simulations allow us to construct statistics on the 
relative difference between the average number of galaxies in the 
North and South Galactic areas of the survey, for each of the three 
samples under consideration. On the observed data, we find that the 
average density mismatch between t he North and the South r egions 
amount to -3.7 %, -0.8 % and -1.9 %. IRoss et al]fe01lll2012l) found 
similar differences in the North/South number densities, and found 
that these differences were consist ent with the mean North/Sout h 
photometric offsets determined by Schlafl v" & Finkbeiner] d201 ll) . 
None of the LOWZ mock simulations provide such a high mis- 
match, while 21 % and only 0.2 % of the mocks are, for the CMASS 
and LOWZ+CMASS samples, at similar or higher levels of North- 
to-South asymmetry. Hence, we conclude that the LOWZ sample 
shows an excessive number of sources in the south part of the sur- 
vey whose effects, however, should be largely mitigated in our anal- 
ysis after removing the dipole in the fraction of area under cosmo- 
logical analysis. 
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where C'f refers to the theoretical angular power spectrum multi- 
pole for our ACDM model, and the superscript ( denotes the par- 
ticular model considered in our two parameter (b, f^ [ ) space. 
This expression is equivalent to the assumption that the covariance 
matrix is dominated by cosmic variance. After taking a flat prior 
for Z^™ 1 in the range f^f e [-300, 300], we compute marginal- 
ized probabilities for both b and ./j^" 1 . We recover practically the 
same central value for the Gaussian bias as in the x 2 analysis 
outlined above. However, the allowed range shrinks considerably 
(b = 1.86 ± 0.09 at 95 % C.L.). On the other hand, the marginalized 
probability of the local non-Gaussian parameter peaks on 59 and is 
contained, at 95 % C.L., in the range /^° L cal 6 [-17, 134]. The width 
of this allowed range increases by less than 10 % when assigning 
the band power centred on multipole £ = 5 a 30 % higher uncer- 
tainty due to systematics (f^f 6 [-20, 143] at 95 % C.L.). We 
next consider including the band power centred on £ = 2, despite 
the concern of it being affected by the dipole subtraction. In this 
case we allow systematics to increase the measured quadrupole er- 
ror by 50 %, and still the allowed range for b does not change, while 
/^° L cal shrinks to f^ 1 e [-2, 77] at 95 % C.L., with the marginalized 
probability of /^° L cal peaking on 38. We conclude that the SDSS DR8 
LOWZ+CMASS galaxy sample is compatible with the Gaussian 
scenario. 

The allowed ranges fo r /^° cal quo t ed her e are in general nar- 
rower than those quoted in Ros s et alj d2013l) (f^ 1 6 [-45, 195] 
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